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Abstract 

An energy-momentum carried by electromagnetic field produced by two point- 
like charged particles is calculated. Integration region considered in the evaluation 
of the bound and emitted quantities produced by all points of world lines up to 
the end points at which particles' trajectories puncture an observation hyperplane 
ip = t. Radiative part of the energy-momentum contains, apart from usual inte- 
grals of Larmor terms, also the sum of work done by Lorentz forces of point-like 
charges acting on one another. Therefore, the combination of wave motions (re- 
tarded Lienard-Wiechert solutions) leads to the interaction between the sources. 



1. Introduction 

We consider a closed system of two point electric charges and their electromagnetic field. A 
charge Ca produces an electromagnetic vector potential A" that satisfies the wave equation 

= -Anj: (1.1) 

together with the Lorentz gauge condition daA^ = 0. The vector is the charge's 
current density which is zero everywhere, except at the particle's position it is infinite. 
For concreteness we imagine that the particles are asymptotically free in the remote past. 

The dynamics of electromagnetic field is governed by Maxwell equations with point- 
like sources. The action of the field of one source on another is described by Lorentz 
force. The evolution of a-th particle is determinated by the relativistic generalization of 
Newton's second law where loss of energy due to radiation is taken into account. 

The dynamics of the entire system is governed by the action 
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Figure 1: The regularization procedure can be performed in two different ways: (i) one when 
Green's functions are used in variational equations of motion; (ii) the other when wave solutions 
are used in Noether conservation laws. 

where /^j, = J2a('^^J■^a■v ~ '^u^a.^I) ■ (a-th point particle carries electric charge and moves 
on a world line (a described by functions z^{Ta), in which Tq is an evolution parameter; 

:= dz^/dTa-) Variation on field variables yields the Maxwell equations. Lienard- 
Wiechert fields are the solutions of Maxwell equations with point-like sources. 

Since the field fa^^u '■= d^Aa^^ — d^Aa^^ generated by a-th source has a singularity on its 
world line, demanding that the total action (1.2) be stationary under a variation 5z^{Ta) 
of the world line does not give sensible motion equations. To make sense of the retarded 
field's action on the particle we should perform the so-called renormalization procedure. 
It involves manipulation of the divergent self-energy of a point charge. As usual, the 
infinite Coulomb-like term is linked with the "bare" mass nia, so that the renormahzed 
mass of particle is considered to be finite. 

The principle of least action (1.2) is invariant under ten infinitesimal transformations 
which constitute the Poincare group. According to Nocthcr's theorem, these symmetry 
properties imply conservation laws, i.e. those quantities that do not change with time. 
In his classical paper [1], Dirac used retarded Lienard-Wicchert solution in the law of 
conservation of the total four-momentum of a composite (one particle plus field, its own 
and external) system. It provides the foundation for his derivation of the radiation- 
reaction force. Lopez and Villarroel [2] substitute the retarded Lienard-Wiechert field in 
the angular momentum conserved quantity which arises from the invariance of the system 
under space rotations and Lorentz transformations. The authors arrive at the angular 
momentum balance equations which is consistent with the Lorentz-Dirac equation. 

To find out Noether quantities G"^ carried by electromagnetic field we integrate the 
Maxwell stress-energy tensor and angular momentum tensor density over a space-like 
three-surface [3, 4, 5, 6]. We obtain terms of two quite different types: (i) bound, G^^, 
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Figure 2: The bound term G^^^ and the radiative term G"^^ constitute Noether quantity 
carried by electromagnetic field. The former diverges while the latter is finite. Bound component 
depends on instant characteristics of charged particles while the radiative one is accumulated 
with time. The form of the bound term heavily depends on choosing of an integration surface 
S while the radiative term does not depend on E. 

which are permanently "attached" to the sources and carried along wdth them; (ii) ra- 
diative, which detach themselves from the charges and lead independent existence 
(see Fig. 2). Within regularization procedure the bound terms are coupled with energy- 
momentum and angular momentum of "bare" sources, so that already renormalizcd char- 
acteristics Gp^J.^ of charged particles are proclaimed to be finite. Noether quantities which 
are properly conserved become: 

= Gp^^i- (1.3) 

Recently [4] a frontal collision of two asymptotically free charges has been considered. 
We have calculated how much electromagnetic field momentum and angular momentum 
flow across hyperplane Et = {y e M4 : if — t}. The crucial issue is that the Maxwell 
energy-momentum tensor density of entire system 

is the sum of individual "one-particle" densities and an "interference" term: 

T^'^ = Tl^^ + Ti;^ + TZ. (1.5) 

An intrigue feature is that the radiative contribution from the combination of the retarded 
Lienard-Wiechert fields 

^^TZ = /s/r2)A + /g/ri)A - 1/4<'^ (/g/i? + /g/ii^) (1.6) 

is then nothing but the sum of work done by Lorentz forces of point-like charges act- 
ing on one another. Therefore, an interference of outgoing electromagnetic waves in an 
observation hyperplane leads to the interaction between the collided sources. (The dif- 
ferentiation of energy-momentum conserved quantity gives the relativistic generalization 
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of Newton's second law [5].) This observation gives us an alternative interpretation for 
the label "int": it stands for "interaction" as well as "interference". 

In this paper we study a closed system of two arbitrarily moving point-like charges 
which are asymptotically free in the remote past. The expressions for work done by 
(retarded) Lorentz forces will be obtained via the rigorous integration of interference 
parts (1.6) of energy and momentum densities (1.5) over three-dimensional hyperplane 



2. Preliminaries 

We choose metric tensor rim, = diag(— 1, 1, 1, 1) for Minkowski space M4. We use Heaviside- 
Lorentz system of units with the velocity of light c—1. Summation over repeated indices 
is understood throughout the paper; Greek indices run from to 3, and Latin indices 
from 1 to 3. The particles' coordinates, velocities etc are labelled a or b. 

We consider an arbitrarily moving particles which are asymptotically free in the remote 
past. Average velocities are not large enough to initiate particle creation and annihilation. 

We suppose that the components of momentum four- vector carried by electromagnetic 
field of particles are [7] 

P:mit) = p[ da,T^\ (2.1) 

where da^ is the vectorial surface element on a observation hyperplane = {y e M4 : 
y^ — t}. Particles' world lines 

t^(t,4(t)) (2.2) 
are meant as local sections of trivial bundle (M4, i, M) where the projection 

i : M4^R 

iy',y')^y' (2.3) 

defines the instant form of dynamics [8] . 

By T'^" we denote the components of the Maxwell energy-momentum density (1.4) 
where field strengths f^" are the sum of the retarded Lienard-Wicchert solutions f^^^ 
and /^2) associated with the first and second particles, respectively. So, the total elec- 
tromagnetic field stress-energy tensor (1.4) becomes the sum (1.5) where the T^'^^ term 
is given by the expression (1.4) where "total" field strengths are replaced by "indi- 
vidual" ones fl^^y The interference term (1.6) describes the combination of the outgoing 
electromagnetic waves. 

The components T^*^ have singularities on particles' trajectories. In equations (2.1) 
capital letter P denotes the principal value of the singular integral, defined by removing 
from an £a-sphere around the a-th particle and then passing to the limit — > 0. 
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Figure 3: The past light cone with vertex at point y G St is punctured by the world lines of 
the 1-st particle and the 2-nd particle at points zi(ii) and -22(^2), respectively. The vector 
is a null vector pointing from Za{ta) = {ta,z1{ta)) to y. 

3. "Interference" coordinate system 

The main goal of the present paper is to compute the interference parts of Poincare 
group conserved quantities carried by radiation. To perform the volume integration an 
appropriate coordinate system for flat space-time is necessary. 

3.1. Local expressions 

The interference terms of energy-momentum and angular momentum at point y G M4 
depend on the state of the charges' motion at the instants ti and ^2 at which their world 
lines intersect the past hght cone (see Fig. 3). Coordinates of an observation point y are 
given by 

y'^^z:{ta)+K: (3.1) 

where is the null vector pointing from Za(ta) G Ca to y. Our next task is to find 
out local expressions for the "light-cone mapping" [9] pictured in Fig.3. We generalise 
coordinate system presented in [4] where a frontal collision is considered. 

The set of curvilinear coordinates contains the "laboratory" time t as well as both the 
"retarded" times ti and ^2- The "laboratory" is a single common parameter defined along 
all the world hues of the system. To find out local expressions for the components of 
null-vectors Ki and K2 we consider an interference of outgoing electromagnetic waves in 
hyperplane (see Fig.4). By this we mean the intersection of spherical fronts S'i(Oi, t—ti) 



5 




Figure 4: The sphere Si{Oi,t — ti) is the intersection of the future hght cone with vertex at 
point zi{ti) G (i and hyperplane S^. The sphere 82(02, t — ^2) is the intersection of and the 
forward light cone of 22(^2) G C2- Intersection S'inS'2 is the circle C{0, h) with radius \OH\ := h. 
It contains an observation point y G (see Fig.3). 



and 5*2(02, t — ^2) pictured in Fig.4. It is the circle C(0, h) centred at point 

{ti-t2){2t-ti-t2) 
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Zl(tl) +Z2(t2) 



+ 



2q^ 



Zl(tl) - Z2(t2) 



(3.2) 



Since = |Z — zi| and IOO2I = |Z — Z2I, the square of the radius h of the circle can 

be expressed in the following alternative ways: 



^ {t-tif -\Z-Zi\ 
= (^-^2)^- IZ-Z2I 



(3.3) 



The characteristics of the circle are obtained from analysis of the triangle O1O2H with 
sides \OiH\ — t — ti, \02H\ —t — t2, and IO1O2I = |zi(^i) — Z2(i2)| Q- 

To define the coordinates of the points of the circle we translate the origin at the 
centre (3.2) of the circle C{0, h) and then rotate space axes till new z-axis be directed 
along three- vector q := zi — Z2 (see Fig. 5). Orthogonal matrix 



to — 




— sin 

cos (pq 

1 




sindn 



1 







(3.4) 



q cos 'dq 



determines the rotation. Finally we obtain coordinate transformation locally written as 

(3.5) 



r = t 

y' = Z\t,h,t2) + h(t,h,t2)u;'j(ti,t2y 



6 



where = (sin 93, cos(/3, 0). Polar angle (p distinguishes the points of circle C{0, h). 

To present the local expressions for the a coordinate system centred on an accelerated 
world line of the a-th particle, we rewrite eqs.(3.5) in a manifestly covariant fashion: 

l/" = <(ia) + f^V(ti,t2)A;:'. (3.6) 

Four components 



= t — ta, kl = h sin (f, kl = hcos(p, k^ 



(3.7) 



satisfy the relations (3.3) and, therefore, constitute null-vector ka- Having rotated it by 
orthogonal matrix ^2 with components Qq^ = Q 



ujij we obtain the vector 



Ka pointing from Za{ta) G Ca to y G (see Fig. 3). The orthogonal matrix uj is given by 
eq.(3.4); it rotates space axes of the laboratory Lorentz frame (see Fig. 5). 
Third component of ka is determined by 



2q 



The characteristics |Z — zi| and |Z — Z2I arc obtained from the analysis of the triangle 
O1O2H with sides \OiH\ = t - ti, \02H\ ^ t - t2 and IO1O2I = |zi(ti) - Z2(t2)| q; 
they are pictured in Fig. (5). 




Figure 5: In "momentarily rotating" Lorentz frame z— axis is directed along three- vector q. 
Circle C{0,h) = Si D S2 lies in Oxy plane; it is centred at the coordinate origin (cf. Fig. 4). 
Polar angle if distinguishes an observation point H G C{0,h). Space parts ki and k2 of null 
vectors ki and ^2 are equal to h sin (fi + h cos ipj + kflc and h sin ipi + h cos ifj + /cf^^, respectively. 



3.2. Global mapping 

To cover the sphere Si{zi{ti),t — ti) where ti is fixed we change the parameter ^2- The 
starting point is the solution t2^*(ii) of algebraic equation 

ti-ff^q{ti,ff) (3.9) 



7 




Figure 6: For a given ti the retarded time t2 increases from tl^^h) to tf^'iti). Minimal value 
^2^*(^i) labels the vertex of forward light cone which is punctured by the world line of the first 
charge at a given point (ti, z\{ti)). The world line of the second charge punctures the future 
light cone of this point at point 

which describes the future hght cone with vertex at (^2*^*, -22(^2^*)) (see Fig. 6). The sphere 
5'2(-22^*, t — ^2^*) touches a given sphere Si{zi{ti),t — ti) at point N (see Fig. 7). If parameter 
^2 increases to ^'^^{ti) being the solution of algebraic equation 

tf"" -ti^q{ti,tf'') (3.10) 

the intersection yS'ifl^'f'^'' contains the only point S. Equation (3.10) looks as the equation 
of backward hght cone of (^2'^", 2^2(^2'^'')), but it defines the future light cone with vertex 
at (ti,z\{ti)) (see Fig. 6). The sphere 5*1 becomes the disjoint union of circles C{0,h) — 
Si n 5*2 if the parameter t2 changes from t2^^{ti) to ^'''^{ti). 

Going along the world line of the first charge we arrive unavoidably at the point t^^*(t) 
being the solution of the algebraic equation 

t-tr^q{tT\t). (3.11) 

The forward light cone of this point touches the world hne of second charge at point 

{t,Z2{t)) (see Fig. 8). Light cones of upper vertices do not intersect the second world 
line at all. Spheres Si{zi{ti),t — ti) determined by ti G [t^^*(t),t] constitute the region 
of hyperplane which requires another parametrization. For a given instant ti from 
this interval the point S (see Fig. 9) is associated with the solution t2{ti) of the following 
equation: 

2i-ii-4 = g(ii,4). (3.12) 
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Figure 7: The sphere S2{02^^ ,t — f2*) is the intersection of the future light cone at {^2*, -^2(^2^*)) 
and T,t. It touches a given sphere S\(P\,t — ti) at point N. The sphere S2{02^'" ^t — t'^''") touches 
S\{zi,t — ti) at point S. If retarded time t2 increases from t'^*{ti) to t'^'"{ti) the sphere Si is 
covered by circles C(0, h) = SiD 82- (A circle S\ fl S2 is pictured in Figs.4,5.) 

The point N in this figure is still connected with the solution t'^2^{ti) of equation (3.9). 

So, we construct the global coordinate system centred on the world line of the first 
particle. It bases on the trivial fibre bundle (2.3). A fibre is a disjoint union of retarded 
spheres centred on the world line of the first particle. A sphere is parametrized by 
the retarded time of the second particle and the polar angle. Locally the coordinate 
transformation is given by equations (3.5). 

In an analogous way we construct the coordinate system centred on the world line 
of the second particle. If e] - oo,t^^*(t)] then ti e [tr*(^2), if "(^2)]; if ^2 e [t2^*(t),t] 
then ti e [tl^^{t2),t'i{t,t2)], G [0,27r[. The ends of intervals are defined by the following 
algebraic equations: 



It is worth noting that the functions ti^*(t2) and tf "(ti) are inverted to each other as well 
as the pair of functions tf ^(t2) and t'^2^{ti) (see Fig.ll). For a fixed observation time t 
the functions t'i{t,t2) and ^2(^)^1) are inverses too. 



j.aav J. 
Il — ^2 

2t - t[ - t2 




q{tT\t2) 

g(tf^^2) 
q{t,tr) 
q{t'i,t2). 



(3.13) 
(3.14) 
(3.15) 
(3.16) 
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Figure 8: The forward light cone of touches the second world line at the 

instant of observation. Future light cones of upper vertices do not intersect it at all. For a given 
ti G [ti^^(t),t] the parameter t2 increases from t2^*(ii) to t2(t,ti). The maximal value t2(t,ti) 
labels the vertex of future light cone which touches the forward light cone of {ti,z\{ti)). The 
minimal value of t2 is the solution t2^*(ti) of equation (3.9). 

4. Electromagnetic fields in terms of "interference" 
coordinates 

Electromagnetic field generated by a— th particle is given by [9, eq.(5.2)] 

ya) 

We use sans-serif symbols for the retarded distance [9, 10] 

ra(y) = -?7a/3(y" - Z^{ta))ul'{ta), (4.2) 

and for the null vector Ka rescaled by a factor r^^: 

K-Uy''-^:{ta)]. (4.3) 

■a 

To rewrite expression (4.1) in terms of "interference" curvilinear coordinates consisting 
of the common evolution parameter t, individual times t\ and t2, and angle variable </?, it 
is advantageous to replace proper time by evolution parameter ta- The components of 
particles' 4-velocities Ua and 4-accelerations a^, a = 1, 2, become [7] 

< = laV^aiQ (4.4) 

< = itiVa-VaK + ^lv^, (4.5) 
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Figure 9: For a given ti G [tT^{t),t] the spiiere 5i(0i,t — ti) is a disjoint union of circles 
C{0,h) = SinS2- Their radius h and centre coordinate Z are determined by t2- The parameter 
t2 increases from t^^^h) (circle 52(0^^*,* - i^^*)) to t'2{t,ti) (circle S2{0'2,t- Q); f e [0,27r]. 



where 4-vectors vj^ = (1, fa(^a)), "^a = (0? Ki'^a)) and factor 7^ := [1 — v^] After some 
algebra, using the relation = K^/ra, we obtain 

.(o) _ Va{ta)Ka^p - Vf3{ta)Ka^a 
JaP — 3 Ca 

' a 

a 

where 

C„ = 7a'+(KaVa), r„ = X° - (K„v„). (4.7) 

Having used differential chart (5.2), one can derive the electromagnetic field (4.6) from 
Lienard-Wiechert potential 

" " r„(y) ^^-'^ 

via the relations /^^^ = aJ^^ - aJ^"^. 
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5. Interference part of electromagnetic field 
four- momentum 



Now, we calculate the interference part of the energy and momentum carried by "two- 
particle" electromagnetic field: 



(5.1) 



An integration hypcrsurface St = G M4 : = t} is a surface of constant t. The surface 
element is given by dao = ^/^dtldt2d^p where is the determinant of metric tensor 
of Minkowski space viewed in curvilinear coordinates (3.5). Differentiation of coordinate 
transformation (3.5) yields differential chart 



d/dy^ 
d/dy^ 
V d/dy^ J 



( 1 



/ 1 








\ 










L1J2I ^^22 








^^33 / 




t-ti 


t - 


-t2 


n 


r2 


/isin (f 


/isin (fi 






r2 


hcosip 


hcos ip 


ri 






|Z - zi| 


|Z 


- Z2 


ri 




r2 



(5.2) 





COS(/? 

sin (/? 
h 



J 



\ ( d/dt \ 
d/dti 

d/dt2 



Its Jacobian gives the determinant of metric tensor mentioned above 

rir2 



-9 = 



<1 



(5.3) 



The volume integration (5.1) can be performed via the coordinate system centred on 
a world line either of the first particle 



tf-{t,) t 



J dti J dt2+ J dti J dt2 ' 



or of the second particle 



dip 



rir2 



(5.4) 



t[(t,t2) 



dU 



dti + 



dU 



dti 



tr*(*2) 



2-K 



dp 



rir2 



(5.5) 



The end points of these integrals arise from the interference pictured in Figs. 4- 9. 
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5.1. Interference part of zeroth component 

In this subsection we trace a series of stages in calculation of the volume integral 



Pint 



^^O^int 



(5.6) 



In Appendix A we perform the computation in detail. 

It is straightforward to substitute the components (4.6) into equation (1.6) to calculate 
the interference part of electromagnetic field stress-energy tensor. We obtain the following 
energy density: 



1 ^dVo 



C2 



+ ro- 



C1C2 



rir2 \dt1dt2r1r2 dtiri{r2f ^^2 (ri)V2 (n)^(r2) 



where function 



of^K Ok Ok 
dtidt2 dtidt2^ 



K 



-Q 



(5.7) 



(5.8) 



does not depend on angle variable at all. 

Taking into account the specific structure of the expression (5.7) which contains the 
partial derivatives we rewrite the integrand ^/—gTP^^. as follows: 



47r rir2^oo 



6162 q 



■-T 



92 



mt 



+ 
+ 
+ 



dtidt2 
d_ 
dt2 



qr\r2j dh 

Cl 



C2 



d ( 



qri{r2Y 
1 \ 



dt2 \qr1r2J 



q{riYr2 dt 
C1C2 d 



1 \qr1r2J 



grir2 

_C2 

_g(ri)2(r2)2 dti V?^i(^2)^ 
1 \ 



d 



Cl 



dh \q{'r'iyr2 



dtidt2 \qr1r2J 



(5.9) 



First of all we should perform the integration over (/? (see integration rules (5.4) and 
(5.5)). The crucial issue is that the integral of the bracketed expression (that which is 
proportional to Tq) over ip vanishes (see Appendix A). Hence the integral of (5.9) over 
the angle variable has the remarkable properties of being the sum of partial derivatives: 



2n 



I 



int 



6162 ^ d 

dtidt2 dti 



Bo 



+ 



d_ 
dt2 



(5.10) 



Here 



1 

2^ 



2tt 



2tt 



dip- 



qr\r2 ' 



27r 



C2 



qr\{r2) 



2' 



1 

2^ 



2-K 



d(p- 



Cl 

q[ri)^r2 



(5.11) 
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where and are given by eqs.(4.7). 

It is natural to integrate tlie expression being tiie time derivative with respect to ti 
according to the rule (5.5). The result is 



t'{t,ti) 



J dti j dt2+ j dti j 



t 



dto 



(5.12) 



dhG2[h,tf''{ti 



dtiG2[ti,tl'\h)]+ / dtiG2[ti,t'^{t,t 



Having apphed the rule (5.4) to the expression of type dGi/dti, we obtain 



t 



t[(t,t2) 



dU 



dti + 



tr*(t2) 



dt2 
t 



dti 



dGijtiM) 
dti 



(5.13) 



dt2Gi[tt^{t2)M- j dt2Gi[ff{t2)M+ j dt2Gi[t[{t,t2),t2 

-oo -oo *2'"(*) 

The double derivative involved in eq.(5.10) can be written in the form either 

'dGo' 



or 



_d_ 

dti 

d_ 



dt2 

dGo 
dti 



(5.14) 
(5.15) 



Now we choose (5.14) and add this term to dGi/dti. 

Therefore, the end points are valuable only in the integration procedure either (5.4) 
or (5.5). The retarded instant, f^^ith), and advanced one, fa^{tb)i {a ^ h) arise naturally 
as the limits of integrals. They label the points S and in which fronts of outgoing 
electromagnetic waves produced by ei and 62 touch each other (see Figs. 7 and 9). Triangle 
O1O2H (see Fig. 5) reduces to the line at these moments. 

An essential feature of integration is that the functions ti*(t2) and t^'^^ti) are inverted 
to each other (see Fig. 11). This cicumstance allows us to change the variables in the "ad- 
vanced" integral involved in eq.(5.13). Further we couple it with the "retarded" integral 
ofeq.(5.12). We obtain 

t 

^1-Vi 



I 



dti 



1 - V, 



Gi — G2 



t2=tr*(tl) 



where 



(5.16) 



(5.17) 
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Scrupulous calculation results the terms of two quite different types: (i) this depends on 
all previous evolution of the 1-st charge 



{5.U 



(ii) those determinated by the state of particles' motion at the observation instant only: 



6162 



1 



ti=t 



ti— >— 00 



2[t 



(5.19) 



(see Appendix E, Table 1, left column, first line). The integral (5.18) over world line ^1 
is then nothing but the zeroth component of the work done by " retarded" Lorentz force 
acting on the first charge. 

It is reasonable that, starting with the retarded Lienard-Wiechert solutions, we obtain 
the retarded direct field due the 2-nd charge on the 1-st one. A surprising feature is 
that wc can arrive at the expression for the advanced direct field within the framework of 
retarded causality. E.g., one can perform change of variables (ti^*(t2)i ^2) ^ {ti^t'i'^^ (ti)) in 
the retarded integral involved in eq.(5.13) and then couple it with the advanced expression 
from eq.(5.12): 



/ 



dti 



Go 



I + V2 



(5.20) 



Having integrated (5.20), we obtain the work done by advanced Lorentz force due to the 
2-nd charge plus functions of momentary positions of particles: 



-/ 



2kl 1 + V2 g[l + ^2] 



tl- 



(5.21) 



(see Appendix E, Table 1, left column, second fine). The matter is that the integral of 
advanced force due to 2-nd charge over worldline Ci is intimately connected with integral 
of the retarded force due to 1-st charge over (2' 



*(t) 



j dt2l2^F',^[ff{t2)M- j dh^^'F^,[h,tf^{h)] 
-00 —00 

-1 + (Vi • V2) , 1 , 1 



(5.22) 



-6162 



?[1 + 1/1] [1 + 1/2] q[l + Vi] q[l + V2] 



t2=t 



t2- 



(see Appendix D, eq.(D.ll)). Therefore, the advanced expression can be replaced by the 
retarded one plus functions of momentary positions of particles (see Appendix E, Table 
2, left column, second fine). 
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Now we consider the last terms in both eq.(5.13) and eq.(5.12). Since the functions 
t[{t,t2) and t2{t,ti) are inverses, the sum of these integrals can be written in the form 
either 

t2=t'^{t,tl) 



dti 



tj<=*{t) 



or 



dU 



l-Vo 



1 + 14 



ii=ti(t,t2) 



(5.23) 



(5.24) 



Both the expressions result the same function of the end points only: 

1 t2^t 



2{t - 12) 



— lim _ , ^ , + 



6162 



,,=,.e*(,) t2-^t2{t-t2) 2[t-V2^\t)] 



(5.25) 



(see Appendix E, Table 1, left column, third line). 

Summing up all the contributions (5.18), (5.19), (5.21) where (5.22) is taken into 
account, and (5.25) we obtain the expression 

t 

dtal-'Fl{t,.tT{Q) (5-26) 



Pint 



6162 



(V2 ■ Vi) + V2 



lim 



6162 V2 



q[t\^\t),t] [1 + Vi] [1 + V2] t2-^t t-t2l + V2 

Now we take the double derivative in the form (5.15) and add it to dG2/dt2. Analogous 
calculations give 



t 

Pit = -E/ dtal-'F,l{t^,tl^\Q) 



(5.27) 



6162 



(vi ■ V2) - 14 ^ e,e2 



ti^t t-til-Vi 



q[tA'\t)] [1-1/1] [I-V2] 

(see Appendix E, Table 1, right column). 

Having compared eq.(5.26) with (5.27) we are sure that the calculations result the 
"immovable core" which describes the action of the fields due to one charge on another, 
and "changeable shell" which expresses the deformation of electromagnetic "clouds" of 
charged particles due to mutual interaction. Only the immovable terms should be taken 
into account in the total energy balance equation. 



5.2. Interference part of space components 

To calculate interference part of electromagnetic field momentum p*^ we have to 
integrate the expression 

^^T^nt = fd)fUi + /(?)/(\), (5.28) 
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over three-dimensional hyperplane ip = t. The electromagnetic field components are 
given in Section 3. 

According to the integration rules (5.4) and (5.5), first of all we perform the angle 
integration. Then integrand (5.28) looks as follows: 



2lT 



6162 



(5.29) 



— ^2 

+ A 

+ Ao 

+ Co 

+ Bo 

+ Vo 



dtidt2 



+ 



+ a 



d^X d^X 



dtidtl 



+ 



dtl 



d^X 
dt\dt2 



k^- 

^dhdt2 ' dtl 



dX 

x{vi+vi) + kX— + ky,:-- 



,dX 



d^X d^X 
^dqdf2 ^ 



d^X 
dtldtl 
d^X 



dt2 

\yi +k%'— + k%' — 
XVl+klV,^, + k2Vi^^ 



■dtl 



dtl 



where 



A = 1/2 [(A;? - k',r - q'] ■ 
CaUigraphic letters A, B, C and V denote the following integrals over (/?: 



(5.30) 



271 



2lT 



A, 



1 



CIC2 



Bl 



C2 



* " 27r ./ "^^^''g(r02(r2)^' ~ L J ''''''' qrMY' 





27r 



2-K 



d^Kl 



Cl 



?(ri)V2' 



I [ 1 

Vl^— d^Kl- 



27r 



(5.31) 



where and are given by eqs.(4.7). Functions BoiCo and "Dq are defined by eqs.(5.11) 
and function ^0 is 



2tt 



Ao 



dip 



C1C2 



2-K J g(ri)2(r2)^ 



(5.32) 



After some algebra one can rewrite the terms which involve Al,Bl, Cl and VI (see the 
1-st and the 2-nd fines of eq.(5.29) as follows: 



_d_ 
dfi 



d_ 

dBi_dCl d^Vj 
dtl dt2 dtidt2 



^b(Cl 



dVl 
dtl 



+ 



dtidt2 



(5.33) 
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Here 



Ai 



9A 



A2 = fc? 



+ 



dX 



(5.34) 



dtidt2 dt2 

Routine scrupulous calculations performed in Appendix B explain that the "non-derivative 
tails" in eq.(5.33) are proportional to three- velocities: 



A.2 



dB\ dC\ 



dti 



dt2 
da 



+ 



dtidt2 



dti dt2 dtidt 



= vi{t2)(c^- 



dVo 
dt2 



(5.35) 



We add them to the part of integrand (5.29) which involve "zeroth" functions Ao^Bq^Cq^ 
and Vq. It is now straightforward (but tedious) matter to rewrite it as the following sum 
of partial derivatives: 



+ 



+ 



+ 



d_ 
d_ 
d_ 



AUBo 



dt2 



d_ 
d_ 

^dr2 



A2 { Co 



dtidt 



0t2 

a(.i + .^)(bo 

[A {v{ + v'2) ^0] 



Al Co 



dVo 
dti 



dtidt2 



+ 



d_ 



X{v{ + vl)(Co 



dtidt2 



dti 
d{v\AiVo) 
dto 



dti 



(5.36) 



(We keep in mind the identity (A. 41)). Recall that Ai, A2 are given by eq.(5.34) and 



Al 



""^^^ dt 



dX_ 
2 



A^ = vik'i 



dX 
dfi 



(5.37) 



Therefore, the integrand (5.29) also becomes the combinations of partial derivatives 
with respect to time variables, namely the sum of the expressions written in the first 
line of eq.(5.33) for both b — 1 and 6 = 2, and eq.(5.36). Now we apply the integration 
procedure developed in previous subsection. 

Each double derivative involved in (5.29) can be integrated according to the rule either 
(5.4) or (5.5). There are five terms of this type in this expression. This circumstance 
implies ten possible ways of integrations. In Appendix E we study two of them in detail 
(see Table 2 and Table 3). 

Firstly we write all the double derivatives in the form (5.14). The integration gives 



int 



+ 



6162 [-l + (Vi.V2)]n^[ir*(0,*] 



(5.38) 



6162 



-\- lim 



q[t\^\t),t] + + 

6162 {r^^[ri^\t2)M-vi{t2)]V2 



q[tl^\t),t] l + Vi 



t2-^t t — t2 



1 - {V2f 



18 



Secondly, we express all the mixed derivatives in the form (5.14). We obtain 

(5.39) 

tl^tt-h 1-(V1)2 

Comparing cq.(5.38) with cq.(5.39), we are sure that the form of functions of momentary 
positions of particles heavily depend on the method of integration. It reinforce our con- 
viction that the changeable " shell" expresses the deformation of electromagnetic " clouds" 
of "bare" charges due to mutual interaction. Thus only the immovable "core" , i.e. sum of 
work done by Lorentz forces of point-like charges acting on one another, possesses relevant 
physical sense. 




ei62 [-l + {^r,■^r,)]nl[t,tr{t)] 



6. Conclusions 

Inspection of the energy-momentum carried by the electromagnetic field of two point-like 
charged particles reveals the essence of renormalization procedure in classical electro- 
dynamics. Volume integration of Maxwell energy-momentum tensor density over three- 
dimensional hyperplane — t gives terms of two quite different types: (i) these depend 
on the state of the particles' motion in the vicinity of the instant of observation t; (ii) 
those depend on all previous time development of the sources. The former involves di- 
verging quantities while the latter contains finite terms only. Structure of the quantities 
which are accumulated with time does not depend on choice of integration three-surface 
while the form of "instant" expression heavily depends on the way of integration. 

"Instant" terms are permanently attached to the charges and are carried along with 
them. By this we mean that a charged particle cannot be separated from its bound 
electromagnetic "cloud" which has its own 4- momentum [11]. This quantity together 
with 4- momentum of "bare" charge constitute the finite 4-momentum of "dressed" charged 
particle. (Note that the electromagnetic "clouds" of sources are deformed due to mutual 
interaction.) All diverging quantities have thus disappeared into the process of energy- 
momentum renormalization. 

The terms which are accumulated with time lead to independent existence. They 
constitute the radiative part of energy-momentum carried by "two-particle" field. It 
consists of the integrals of individual Larmor relativistic rates over corresponding world 
lines and the work done by Lorentz forces of point-like charges acting on one another. 

The situation considered here, in which the radiation is propagating outward, breaks 
the time-reversal invariance of Maxwell's theory. Choosing the retarded solution of wave 
equation (1.1) as the physically-relevant solution, we adopt a specific time direction, 
when an interference of outgoing electromagnetic waves leads to the interaction between 



19 



+ 00 




Figure 10: To restore time-reversal invariance we locate the observation hyperplane y° = t in 
the distant future. We suppose that particles are asymptotically free. 



the sources. The interference are pictured in a fixed observation hyperplane = G 
M4 : y° = t}. To restore time-reversal invariance we take the limit t +00 and suppose 
that particles are asymptotically free in the distant future. The relation (D.12) takes the 
form 



+ 00 +OD 



(6.1) 



The work done by retarded Lorentz force of 6-th charge over entire world line of a-th one 
is equal to the work done by advanced Lorentz force of a-th particle acting on 6-th charge 
backward in time! The sum of " retarded" works involved in the total energy-momentum 
of our closed (two particles plus field) system 



+00 



+CX) 



j dt,-f^'F^,[tiA'\ti)\ + j dt2i^^F^^[tr\t2)M, 

—00 —00 

may be replaced by the linear superposition 

" +00 

—00 
+00 



(6.2) 



+ J dt2l^' {F^2[fi\t^), t^] + ^^^[^(^2), ^2]) 



(6.3) 



which restore time-reversal invariance. Indeed, the retarded Lorentz force F^[ta, *r*(^a)] 
becomes the advanced one F^J(ta, t'^'"{ta)\ (and vice versa) if the time direction is reversed. 
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But the retarded causality is still not violated. We consider the interference of outgoing 
waves at distant future instead of a picture in which the radiation is propagated inward. 

The situation looks as that described by Wheeler and Feynman [12] where the absorber 
theory of radiation is elaborated. The basic assumption is that the fields which act on a 
given particle are represented by one-half the retarded plus one-half the advanced Licnard- 
Wiechert solutions of wave equations. To disappear "incoming" radiation, the authors 
introduce a perfect absorber which cancels the (acausal) advanced part of the fields acting 
on a given particle and doubles the retarded one. 

Our emphasis is on rigorous calculations and exact solutions based on standard electro- 
dynamics. It allows us to substitute the phenomenon of interference of outgoing electro- 
magnetic waves for acausal mechanism of perfect absorbtion in time- symmetric action-at-a 
distance electrodynamics. The interference of outgoing electromagnetic waves (retarded 
Lienard-Wiechert solutions) ensures the action of the field of one source on another (mu- 
tual interaction). 
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A. Angle integration of 7]°^ 

In this Appendix we perform the integration over (f of "double zeroth" component (5.9) of 
the Maxwell energy-momentum tensor density. Angle-dependent terms involved in energy 
density have the form 




(A.1) 





















where the retarded distances are 



Ta = r° — sin — r^ cos (p. 



(A.2) 



The other scalars we use are: 



Ca = -c° + c^sin(/? + c^cos(^ 



(A.3) 



Here 



r'a - k^-iv^.nX, rl^ 




^1 = ^2jvih 



(A.4) 
(A.5) 
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It is convenient to introduce three-dimensional manifold Q with space-favouring metric 
Qap — diag(— 1, 1, 1). For Q we put the tangent bundle TQ being the disjoint union of all 
tangent spaces T^Q- A tangent vector with foot point a e Q is simply a pair (a, r) with 

r = r"e„ e M^ 

where := d/dx'^^a = 0,1,2, is the standard basis of R^. We define also cotangent 
bundle T*Q being the disjoint union of T*Q. An one-form with foot point a e Q is a pair 
(a, r) with 

r = rpef^ E M^ 

where e^, /3 = 0, 1, 2, constitute dual basis e'^(ea) = 5^. We shall use ga/3 = diag(— 1, 1, 1) 
and its inverse g"^ — diag(— 1, 1, 1) to lower and raise indices, respectively. 

For each differential manifold one can define the canonical pairing 

(,) : T*QxqTQ^R (A.6) 
(ri,r2) ^ ri,^r] 

where both one-form f i and vector r2 are of the same foot point. We introduce also the 
scalar product 

(•) : TQxqTQ^R (A.7) 
(ri • rs) ^ gapTiT^ 

which is connected with canonical pairing by the operation of raising indices. And finally, 
we shall need the norm of vector r e TQ 



-9apr^rf^ (A.8) 



So, a-th retarded distance becomes the scalar product of the vector with components 
(A. 4) and the null-vector := (l,sin</?, cos</?) taken with opposite sign. 

To go further we express a term of type (b • n^)/rir2 as follows 

(A.9) 



rir2 



-A-C\: 




\cos{ip + l3i) ^ -B + C\r2\ 


cos{(p + (32) 


'1 1 






k2 





where scalar b denotes the product (b • n,^). The a-th phase is determined by the 
relations 

cos/3„ = r^/|r„|, sin/3„ = r^/|r„|, (A.IO) 

where 

= h^vl - (v„ . n,y. (A.ll) 
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The coefficients A, B and C are the solutions of the following system of algebraic equations 



(C A B) 



( Lq L\ L2 ^ 

r2fl r2,i r2,2 
\ ri,o ri,i ri,2 / 



(A.12) 



where Lq, = QapL^ and is '^-th component of the vector 

^ eo ei 62 

L= r2,o r2,i r2,2 
V ^^1,0 ri,2 



(A.13) 



This can be written more compactly 



by use of the Ricci symbol in three dimensions: 

+1 when q;/37 is an even permutation of 0, 1, 2 
—1 when q;/?7 is an odd permutation of 0, 1, 2 
otherwise 



-a/37 



(A.14) 



In solving the problem (A. 9) one is soon led into rather complex expression. Great 
simplification arise, however, when one uses the binary operation of vector product which 
is defined as follows: 



[x] : r*QxQT*^ 



TQ 



(A.15) 



So, the determinant D of 3 x 3 matrix in eq.(A.12) becomes the square of vector product 
of one- forms fi and f2 given by eqs.(A.4): 

D = (L,L) = (L)2 (A.16) 

= [f2 X 

= -(f2,r2)(fi,ri) + (f2,ri)2 
= -(r2 • r2)(ri • ri) + (r2 • ri)2 

Having solved the system of linear equations (A.12) we obtain 

(b-[fixL]) 



D 



(A.17) 



B 



--^ |(b ■ r2)(ri ■ ri) - (b ■ ri)(ri • r2; 
(b • [f2 X L]) 



D 



i (b • r2)(r2 • ri) - (b • ri)(r2 • V2) 



C 



(b-L) 

D 



(A. 18) 



(A.19) 
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The expression of type (A. 9) can be integrated over (p via the relations 

' . = "'^^ =0 (A.20) 



1 — asin(/7 — a2 ' 1 — a sine/? 

< a < 1. 

Having integrated (A. 9) we obtain 

/ = -2n —— + —— . (A.21) 

Having considered the simplest case of integral Vq (see eqs.(5.31)), we put the one-form 
b = (-1,0,0). We obtain 

I,„ = -1 + , (A.22) 

5 Vl|ri|| ||r2||/ 



where 



A, = f!!!^ (A.23) 
r^(ri • ri) - r?(r2 ■ ri) 



5n 



Co 



D 

D 

rl{v2 ■ ri) -rg(r2 • V2) 
D 



To calculate Bq we rewrite the integrand as follows 

C2 _ C2 1 

ri(r2)^ rir2r2 



-^2 - C'2|ri| cos((^ + (3i) ^ -B2 + C2|r2| cos((^ + /92)^ 1 





C'2\ 




COs((/9 + j3i) 


' 1 


|ri| sin(v9 + /^i) 
C2 r2 cos(v? + 132) 


[^2° - 1 


r2| 


sin((^ + /32)]^ 



ri r2 J r2 

+ 



-B2 + C2 


r2 


COs((/9 + /?2) 


'2 


r2 


sin(93 + (32) 



where 



(C2 ■ [fi X L]) (C2 ■ [f2 X L]) (C2 • L) 

^2 - ^ , -D2 — , O2 - — , {A.Zb) 

A'^ = -A2Ao + C2Co||ri||2, S^ = -^2So + C2Co(r2-ri), 

= -A2C0-C2A0 (A.26) 
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Another relevant integration rules are 

/ ^^^71 ■ ^i = n / ^'^Tl = ° ^"^^ 

Jo [1 — asmip)^ (1 — a^)'^/^ Jq (1 — asin(^)^ 

< a < 1. 



Bo = -- ( ^ + ^ + ) . (A.28) 



Combining these results together with the relations (A. 20) for integral of expression (A. 24) 
over ip gives 

? Vl|i"i|| ||r2|| ||r2| 

If one interchanges the indices "first" and "second" in the above expression (A. 24), 
they obtain 

1 / A\ B\ Air? \ , , , 

where 



(ci • [ri X L]) (ci • [f2 X L]) (ci ■L) _ 

A = ^ , = ^ , = — — , (A.dUj 

= -5iAo + CiCo(r2Ti), Si = -SiSo + CiCo||r2||', 
C( = -SiCo-C^i^o- (A.31) 

We now turn to the calculation of ^o- Transformation of the integrand scales as 
{Tir2)~'^ proceeds with the help of eq.(A.9), using identity 

^^^^ (A.32) 



The calculation is straightforward, although it involves a fair amount of algebra. Finally 
we obtain 

ci C2 _ /i + i"o|ri|cos((/? + /3i) ^ Ji - Jolral cos((/? + /32) ca QQ^ 



rir2 rir2 {nY (ra)^ 

_^ -^12 - Ci2|ri| cos((/? + /3i) ^ -B12 + Ci2|r2| cos((/? + P2) 
ri r2 



where 



h 


= A1A2-C1C2II 


rif, /o 


= A1C2 + A2C1, 


Jl 


= B1B2 — C1C2I 


r2|^, Jo 


= B1C2 + B2C1, 


A12 


= [A1B2 + B1A2 


-2CiC2(ri 


■ rs)] Ao - JoCo||ri||2 - /oCo(ri ■ rs) 


B12 


= [AiB2 + Ba2 


- 2CiC2(ri 


• ra)] Bo - JoCo(ri • r2) - /oCo||r2p 


Ci2 


= [A1B2 + B1A2 


- 2CiC2(ri 


■ r2)] Co + IqBo + JoAo. 



(A.34) 



(A.35) 

Using integration rules (A. 20) and (A. 27), we perform the integration of (A. 33) over the 
angle variable (p: 
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All the coefficients involved in resulting expressions (A. 22), (A. 28), (A. 29), and (A. 36) 
should be rewritten in terms of three-dimensional vectors which denote particles' positions, 
velocities and accelerations. Substituting components (A.4) and (A. 5) into expressions 
(A. 23), (A. 25) and (A. 30) returns the root coefficients Ai, and Cj, i = 0, 1, 2: 

^0 = _(KxviM".xl]), (A.37) 

(K X V2] • [n, X 1]) 
rSn — 



where 



Co 



A 

^n, ■ [V2 X vi]) 



A = [n, X If - /i'(n,[v2 x Vl])^ 1 = r°Vi - r°V2; (A.38) 

A, = A-^ |c? (K X vi] ■ [n, X I]) - r? ([n, x vi] ■ [n, x 1]) 

+ /i' (n, • [vi X vi]) (n, • [vi x V2])| , (A.39) 

B, = -A-i |c? ([n, X V2] • [n, x 1]) - ([n, x Vi] • [n, x 1]) 
+ /i^ . X V2]) (n^ • [vi X V2])| , 

Ci = A-^|-c?(n,-[v2XVi]) + (n,-[vixl])|, 



(A.40) 



A2 = A-i |c° ([n, X vi] • [n, x 1]) - r° ([n, x V2] • [n, x 1]) 

+ /i^ (rig • [v2 X vi]) (iiq • [vi X V2]) I , 

B2 = -A-^|c°([n, xv2]-[n,xl])-r2°(Kx V2]-Kxl]) 

+ (iiq • [v2 X V2]) (iiq • [vi X V2]) I , 

C2 = A-i|-c°(n,.[v2XVi]) + (n,-[v2xl])|. 

A complex calculation performed with the help of software system " Maple 8" confirms 
the key identity 

^„_|._|. + ||^.0. (A.41) 

OTi dt2 OT1OT2 

It allows us to rewrite the integral of "double zeroth" component of the Maxwell energy- 
momentum tensor density over (/? as the sum (5.10) of partial derivatives in time variables. 



26 



It is worth noting that all the coefficients (A. 37)- (A. 40) and, therefore, expressions 
(A.22), (A. 28), (A.29) and (A.36) depend on h'^, i.e. on the square of the radius of the 
circle C{0,h) — Sid S2 (see Figs. 4, 5). One can express functions ^O: Co and T>q in 
form of expansions in powers of h"^. (To simplify the calculations as much as possible we 
can rewrite the integrands of (A.l) as expansions in power h and then integrate over ip.) 
Putting /i^ — > we tend to convex neighbourhood of the end points, either S or N (see 
Figs. 7, 9). The identity (A.41) is also valid in the immediate vicinity of the end points. 
(Differentiating functions Bo, Co and Vq in time variables we must keep in mind that 
dh'^/dta does not vanish even if h"^ 0.) 

B. Angle integration of T^^ 

To express the integral of 7]^^ over </? as a combination of partial derivatives in time 
variables we have to calculate the following "tails": 

(see eq.(5.33). By calligraphic letters we denote the integrals over angle variable: 



2tt 2tt 

1 f , ci If, Ki 



CI = — d^Kl—^, vi = — 



2tt J g(ri)^r2' 2tt J qrir2 





where 



Kl^ klug + hujl sin (p + hu;2 cos (fi. (B.3) 
The integration can be handled via the relation (A. 9). The simplest term T>1 becomes: 

Here 



Q F2I 



4 = klniAo - ^ { [4 - <(n,V2)] (ri • 

- [^i-n;(n,v0](ri-r2)}, (B.5) 
Bi^klniBo + ^ I [^^ - n;(n,V2)] (r2 • ri) 

- bi-<(nqvi)] (r2T2)| 

where Aq, Bq and A are defined by eqs.(A.37) and (A. 38). But we find out the expressions 
(B.l) in another way. 
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To simplify the calculations as much as possible we express the integrands of eqs.(B.2) 
in form of expansions in powers of h. Thanks to exponential operator 



Y 



exp 



d 



- XI sin (p + UJj2 COS ip) — 



(B.6) 



we remove harmonic functions from denominators and then integrate over if. In fact, we 
deal with the flow of the vector field in between the square brackets of eq.(B.6). It maps 
an open neighbourhood of end points either S or N to an open vicinity of another point 
of integral curve of this vector field [13]. It is sufficient to compute "tails" (B.l) at the 
end points where h'^ — (see Figs. 7 and 9). 

At these end points the term Al is as follows: 



4 



27r 



C1C2 







?(ri)2(r2)^ 



(B.7) 



^2=0 



„0„0 



Since derivatives dh'^/dta, a = 1, 2, do not vanish whenever h"^ 
and Bl up to the first order of this small parameter: 



+ 



kln\BQ 
2gr0(r^)2 



-c^ ( 2- 



ni{nqV2) v{ - n^(ngVi) 



+ 



0, we should expand CI 
(B.8) 



(B.9) 



+ 



^0 
'2 



'1 



Symbols Bq and Co denote the expansions of corresponding integrals (5.31) in powers of 



Bo - 

+ 
+ 

Co = 



9^1(^2)^ 



2grO(rO)2 
1 

-c° 



Kva]^ , ^([ngVi][ngV2]) , [n^vi] 



+ 2 



o„o 



+ 



1' 2 



o (Kv2][ngV2]) , ([ngV2][ngVi]) 

n ~r n 



(r?)2 



(B.IO) 



fB.lll 
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+ 
+ 



2?(r?)V0 



2g(r?)V0 



[ngV2]2 + o ([ngVi][ngV2]) ^ ^ [ngVi]^ 



([ngVi][ngV2]) , ^ ([ngVi][ngVi]) 

n ~r ^ n 



».0„0 
' 1'2 



(r?)2 



The last expansion we shall need is 



+ 



+ 



+ 



+ 



2qr^r2 



'1 



[^-n;(ngV2)][ngV2] 



2{v\ - n^(nqV2)]([ngVi][ngV2]) + {v{ - n^(nqVi)][n<^V2]^ 



r0(r0)2 

2{v\ - n^(ngVi)]([ngVi][ngV2]) + 



,(nqV2)][ngVi]^ 



+ 3 



(,0)3 



By Do we denote the following expansion: 



qr\r\ 



+ 
+ 
+ 
+ 



j [n,vi]2 ^ ([ngVi][ngV2]) ^ [ngV2]2 



0^0 



1' 2 



„0\2 



0^0 



1' 2 



(r?)3r0 



2([n,Vi][ngV2])=^ + [ngVi]^[n,V2]' 
(r?)2(r0)2 

„Kv2]-([nyVi][ii,,V2]) „[nyV2]' 



r?(r2°)^ 



+ 3^ 



0\4 



+ 0{h^). 



(B.12) 



(B.13) 



Our final task will be to compute expression (B.l). When we differentiate functions 
Bl, CI and VI we must keep in mind that derivatives of with respect to ta do not vanish 
even if /i^ — > 0. With a degree of accuracy sufficient for our purposes we obtain 



dti 



dt2 



+ 



dtidt2 



dh dt2 dtidt 



dt2 



(B.14) 



Substituting these relations into cqs.(5.33) returns the integral (5.29) of interference part 
of the momentum density T^^^ over (p as the combination of partial derivatives in time 
variables. 
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C. Direct particle fields and Lorenz forces 



In classical electrodynamics the four- dimensional delta function of the square of the inter- 
val between points A and B is Green's function of the wave operator. The delta function 
ensures that the typical points A and B on the worldlines of point-like charges a and b 
interact if and only if they are connectible by a null ray. The interaction is described by 
Lorentz force, i.e. there is no self action. 




Ca Cb 



Figure 11: Points A G Ca and B E (b are connectible by a null ray. They are defined by the 
pair of instants either {ta,tl^^{ta)) or {t^'^'"{tb),ti,)- Functions tl^*{ta) and t^''''"{tb) axe inverses. 

The particle a is acted on by the particle b via Lorentz force F^" = GaF^-^pU^ where 
F(^)fj is direct particle field [14]. By this we mean electromagnetic field generated by 
6-th particle at point where a-th particle is located. It immediately implies /i = in 
expressions (4.6) for the components of electromagnetic fields. Indeed, h is the radius of 
the circle fl S*;,, i.e. of the intersection of spherical fronts of outgoing electromagnetic 
waves generated by charges (see Figs. 4, 5). If we consider the direct particle field, the 
sphere Sa reduces to the point where a-th particle is placed. 

To evaluate the retarded field of the 2-nd particle at point Zi{ti) G Ci we put = 
and kl = t2^*(ti)] in (4.6). It implies 

= ii'^ = gn;, r2 = q[l-V2], C2 = j^' + qV2 (C.l) 

in the expression for f^^. (It is obvious, that fj^^ vanishes.) All the quantities are 
evaluated at the moments either ti or i2^*(ii), V2 :— (n^ • V2) and V2 '■— (n^ • V2). 

To find out the advanced field of the 1-st particle at point ,22(^2) £ C2, we put k2 — 
and k1 — —q[t1'^'"{t2),t2] in f^^ given by eq.(4.6). It means 

K',=-q, Kl = -qni^, n = -g [1 - K] , c, = - qV, (C.2) 

where Vi :— (n^ • Vi) and Vi :— (n^ • Vi). 
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In general, to obtain the retarded/advanced field generated by a-th particle at point 
where &-th particle is located, one should substitute the quantities 

Kl^eq, Ki^ = {-irqni, r„ = eg [1 - (-1) W„] , c„ = 7„-' + (-1)'^?^ (C.3) 

in eq.(4.6). Parameter e is equal to +1 for retarded fields and —1 for advanced ones. 
Putting eqs.(C.3) in (4.6) we arrive at the following expressions: 



-^In +1"-'-) Tl^a + e 9 

g2 (1 - " g (1 - £aK)' g (1 - ^aK)' 



\a a. q a q 



Q (1 - £aK) g (1 - £aK) 



(C.4) 



where parameters 

Sa = (-l)"e. (C.5) 
The components of Lorentz force a-th charge acting on 6-th one are written as follows: 



— e- 



SgVh - (Vg • Vfe; ^,-2 



q{l-eaVa)- 



+ e 
+ 



-e,F^\e)-re,F^f{e)vi 

-&b&a ^£an\ 

g(i-£aK)' 

1 - £aH " 



1 - SaVa 



la + (-1)'^ 



a -l + (Va-n) ^> 



2^ 



g(l-£aK) 



+ e 



g(i-£aK) 



(C.6) 



(C.7) 



All the quantities labelled by a are referred to the instant ^^(tb) while those supplemented 
with index h are evaluated at t},. 



D. Difference of work done by "advanced" and retarded Lorenz forces 

The retarded, and "advanced" , t^'*'(ta), instants arise naturally within the integra- 

tion procedure developed in Section 5 as the end points of "inner" integrals (see eqs.(5.4) 
and (5.5)). Typical points A (on the worldline of charge a) and B (on the worldline of 
charge h) interact if the line connecting them is a null ray. It seems, that the interaction 
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can be both forward {B to A) and backward {A to B) in time (see Fig. 11). And yet the 
retarded causahty is not violated. Indeed, we consider the interference of outgoing waves 
present at the observation time t. Both the retarded and "advanced" moments are before 
t. 

In this subsection we compare the work done by retarded Lorentz force due to charge 
b on charge a 



(D.l) 



and the work done by "advanced" response of charge a on charge b 



(D.2) 



The following identity generalises the derivatives of eqs.(3.9), (3.10), (3.14) and (3.15): 

dtb 1 - EaVa 

Here 

Sa = i-lTe; (D.4) 

parameter e is equal to +1 for retarded instants and —1 for advanced ones. With the help 
of eq.(D.3) we obtain the following chain of identities: 



dq[tl{h),h] 
dtb 

dnl[tl{tb),tb] 
dtb 

djUq ■ Vb) 

dtb 



(-1) 

(-1)" 
Q 

i-iy 



1 - SaVa 



i ^ l-eaVb iVb-Va 

vl - v„- — - n' 



(D.5) 
(D.6) 



-lb - [-1 + (vi • V2)] — — - 

J- ^n.y n. 



+ SaiVb - Va, 



1 - £aK 



(D.7) 



d(n,-v„) _ (-1)" 



dtb 



q 



7a '^-^ - 1 + (Vi ■ V2) + SaiVb - K) 



1 - SaVa 



' 0,-1 -f r 

1 - e„Vn 



(D.8) 



To compare (D.l) and (D.2) we change the variables [ta'^^{tb),tb\ >—>■ [ta,tl^*{ta)] in 
"advanced" integral: 



t "-6 v"-; 

I dta%'Fi:^[t,,ti^\t,)]- I dtb%'F^,[t:'^{tb),tb] 



(D.9) 



-00 
t 



= J dta 



tb=tr\ta) 
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Using identities (D.6)-(D.8) in the integrand of eq.(D.9), we derive that it is the total 
time derivative. In other words, the difference of "retarded" work (D.l) and "advanced" 
one (D.2) is the integral being a function of the end points only. 




Figure 12: Difference of work done by retarded Lorentz force due to charge 2 on charge 1 and 

the work done by advanced response of charge 1 is given by eqs.(D.lO). All the quantities in the 
right-hand side of this equation which are labelled by 1 are referred to the instant of observation 
while those supplemented with index 2 are evaluated at f^^it). 



11 = , -6162 



ti=t 



(D.IO) 



H = I 



-6162 




i-i + (v,v2)j»; ^ 

q[l - \4J[1 - I/2J g[l - Vi\ "~ q[l - V2\ 



tY'''{t) 



—00 

1 + (V1V2) 



(D.ll) 



// = , 6162 



11 = I 



6162 



g[l + 1/1] [1 + 1/2] g[l + 14] g[l + 14: 

" [-1 + (viV2)]n; _ v\ _ v\ 



- t2^— 00 



t2^— 00 
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Figure 13: Difference of work done by retarded Lorentz force due to charge 1 on charge 2 and 
the work done by advanced response of charge 2 is defined by eqs.(D.ll). Ah the quantities 
in the right-hand side of this equation which are labelled by 2 are referred to the instant of 
observation while those supplemented with index 1 are evaluated at V^*{t). 

It is convenient to rewrite the results (D.IO) and (D.ll) in a manifestly covariant 
fashion: 



I 



:-i)%e2 



q{nq ■ ui) iuq ■ U2) q{nq ■ ui) q{nq ■ U2) 



ta=t 



:d.i2) 



Symbols u'^,a — 1,2, denotes the (normalized) four-velocity vector {'^a^ ,%^v]j). If the 
2-nd particle moves in the retarded field of the 1-st one while the 1-st particle moves in 
the "advanced" field of the 2-nd one, then — (l,n*). Four-products of this null vector 
with four velocities are as follows: 



{Uq ■ Ua) = -7a ^ [1 " K] ■ 



(D.13) 



If one interchanges the words "first particle" and "second particle" in the above sentences, 
= (— 1, n*) and we have 

K 7a Ml + K]. (D.14) 

in eq.(D.12). 



E. Time integration of Tj 



int 



In this paper we integrate the interference part of energy-momentum tensor density of 
two point-like charged particles over three-dimensional hyperplane = {y e M4 : y° = 
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t}. An integration hypersurface is a surface of constant value of the obsevation time 
parameter. Besides t, the set of curvihnear coordinates includes the "individual" retarded 
times ti and t2, associated with the particles' worldlines, and the angle variable (f. The 
integration over (/? is performed in Appendix A and Appendix B. The crucial issue is 
that the resulting expressions are the sum of partial derivatives in individual times (see 
eqs.(5.10) and combination of (5.33) and (5.36)). It allows us to perform the integration 
over one of the time parameters, either ti or t2- "Retarded" shifts in arguments of 
particles' individual characteristics such us coordinates, velocities etc. appear on this 
stage as well as "advanced" ones. 

The first double integral involved in the rules either (5.4) or (5.5) defines the integration 
over "causal" region which is pictured in Figs. 6 and 7, while the second one deals with 
"acausal" region (see Figs. 8 and 9). The integration of "causal" type can be handled via 
the relations (D.3)-(D.8). Their counterparts for "acausal" region look as follows: 



dt[{t,t2) 
dt2 

dq[t[{t,t2),t2] 
dt2 

dni[t[{t,t2),t2\ 
dt2 

djnq ■ vi) 
dto 



djuq ■ va) 
dU 



I-V2 
'l + Vi 

V1 + V2 
" 1 + 14 



1 
q 



,l-V2 ,V, + V2 

+ n 



7i 



,1-1^2 



1 + Vi 

V2 



1 + (vi ■ V2)] - 



(E.l) 
(E.2) 
(E.3) 



1^1 + ^2 



(E.4) 



72"' - [1 + (vi • V2)] 



1 + T4 



{Vi + V2) 



Vo 



1 + ^1 



+ V2 



(E.5) 



where V; := (n^ ■ v^) and K 

The way of integration where all the mixed derivatives are written as d / dti[dGQ/ dt2\ 
results the expressions placed in the left columnes of Tables. We apply the rules (5.16), 
(5.20) and (5.24) for the 1-st, 2-nd and 3-rd line, respectively. 

If one changes the order of differentiations they obtain the expressions in the right 
columnes of Tables. For the 1-st and 2-nd fines time integration rules are as follows: 



dU 



1-^2 
1-1/1 



tl (*2) 



a), / dt2 



(b). (E.6) 



ti=tr*(*2) 



Acausal region is integrated according to the rule (5.24) (3-rd line of the right column). 
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Table 1. Integral Jo'' V^T^nt has the form of dGi/dti + OGs/Sis + d'^Go/dtidt2. In- 
tegration over time results the expressions in the left column (if mixed derivative is coupled 
with dGi/dti) or in the right column (if d'^Go/dtidt2 is added to dG2/dt2)- Integration over 
"acausal" region gives the functions of the end points only (see third line). 



d \dGo] 
'dtl [ dt2 J 


d pco 

9^2 L C'^l 




(5.16) 
+6162 


" 1 1 + K2 1 


^*(^i)] 
- ti=t 

. ti— >— 00 


(E.6a) 
+6162 


dt272-'n[^r 

1 1 + 1 1 


{t2),t2] 
t2^—00 


_2A;^ 1 - V2 q[l - V2\ 


2A;i^l-l/i ' q[l-V^\\ 


(5.20) 
+6162 


1 I-V2 1 1 


ti— >— 00 


(E.6b 
+6162 


)-f_^dt2j2'FUtr{t 
' 1 1 


2),t2] 
t2=t 

t2— »— 00 


2k'il + V2 ' q[l + V2\_ 


2k"^l+ Vi q[l + Vi\ 


(5.24) - If 


(5.24) 1^ 


t2=t 



Table 2. Integral Jq^ y/—gT^nt becomes the combination of partial derivatives in time 
variables. Structure of this Table is analogous to the structure of Table 1. 



d \dA] 
'dtl \_dt2. 




d \dA' 




+6162 




] 

■1 tl=t 

. tl— >— 00 


+6162 


r}ydt2^2'Fi2[tnt2), 

ni + v\ vi 


t2] 

t2-^t^''*{t) 
t2— »— 00 


[2kl[l - V2] q[l - V2] 


2kAl-V,] q[l-V^\\ 


+6162 






+6162 


t^dt2l2'F{2[t\^\t2).t2 


] 

1 t2=t 
t2— > — 00 


2k'i[l + V2] q[l + V2]_ 


2k\\\ + V^i] g[l + 


'''' 2e,[l - V2] 


t2^t 


nl + vl 
2k',[l + V,] 


t2=t 



Taking into account the relationship (D.12) between work of the "advanced" Lorentz 
force and the work of the "retarded" one we remove all the "advanced" integrals from 
these Tables. The final expressions are written in Table 3. 
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Table 3. The expressions which are placed above double line concern with integration of 
energy density Tj™ while ones below double line result from integration of T^^j. 



d 

mi 



OA 



m 



2 J 



+6162 



J_l±ii _ 1 
2A;°l-^2 q[l-V2\ 



ti=t 



ti— >— 00 



1 1-V2 1 



-6162 



2kl 1 + V2 q[l + V,] 



-1 + (V1V2) 



-1 t2^t 



6162 

2^2 



t2^t 



+6162 



tx=t 



ti— >— 00 



+6162 



[-l + (viV2)]n; 



t2— > — 00 



2k\{l - V^] 



d 

W2 



OA 

mi 



+6162 



1 1 + Vi 



2k1 1 - q[l - ¥2] 

-1 + (V1V2) "' 
q[l-V,\[l-V2\ 



-lLdt2l2'FUtr(t2)M 



+6162 



1 l-Vi _ 1 
2A;? 1 + V, q[l + V,\ 



t2=t 



t2— »— 00 



6162 

2kl 



t2=t 



-lLdtiiT'F^i[ti,tr(h)] 



+6162 



nl + vl 



2k'l[l - Fi] g[l - 1^2j 

[-l + (viV2)]n; -*^^* 



-lLdt2l2'Fl2[t?\t2)M 



+6162 



2A;i'[l + V,] q[l + Vi] 



t2=t 



t2— » — 00 



K + y 



6162 



2k\[l + 1^1] 
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